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Existing non-perturbative computations of the running coupling of quenched QCD in the Schrodinger functional 
scheme are extended to scales fj, lying much deeper in the low-energy regime. We are able to reach /x" 1 ~ 0.9 fm, 
where a significant deviation from its perturbative evolution is observed. 



1. MOTIVATION 



2. SF SETUP 



The Schrodinger functional (SF) of QCD 
provides a suitable tool to compute the run- 
ning of the strong coupling a s and to 
solve other renormalization problems ^|JJ] non- 
perturbatively by means of numerical simula- 
tions. The basic idea to cover the scales involved 
- often differing by many orders of magnitude - is 
a recursive finite-size scaling technique |5j to re- 
late the low-energy sector with the scaling regime 
at high energies, where the evolution of a s follows 
the perturbative renormalization group. 

A particular property of the coupling in the SF 
scheme, a$F, which however should not be re- 
garded as universal for QCD couplings, is that its 
non-perturbative running is quite accurately de- 
scribed by perturbation theory (PT) down to sur- 
prisingly low energies. Therefore, it appears in- 
teresting to investigate asF(n) in quenched QCD 
at much lower scales \x = 1/L (L the linear size 
of the finite system) than considered so far jl],[| 
in order to see a deviation from PT. 

As an aside we would like to remark that dur- 
ing the year 2000 the present study at the same 
time was designed to test the APEmille comput- 
ers installed at DESY Zeuthen || and to gain 
experiences with code adaptation and optimiza- 
tion for this architecture in the context of a first, 
clear-cut physics project. 



*Based on a poster presented by J.H. at the conference 
LATTICE '01, August 19 - 24, 2001, in Berlin, Germany. 
T Supported by the EU under HPRN-CT-2000-00145. 



The (lattice regularized) SF is given in terms of 
the effective action (free energy) T of QCD satis- 
fying Dirichlct boundary conditions in Euclidean 
time 0. A renormalized coupling, g| F = g 2 , 
is introduced as the response to an infinitesi- 
mal variation of a specific 1-parameter (r/) fam- 
ily of prescribed constant abelian boundary fields. 
Taking into account the perturbative series T = 
To 9 2 + Ti + T 2 gl + ■ ■ ■, we define g 2 by 

&£ = dTo _J_ 
9v v=0 drj 7j=Q g 2 {L) ' 

~g genuinely depends only on one renormaliza- 
tion scale, L = l//i, and the quantity that has 
been devised to map out the scale evolution of ~g 
is the step scaling function (SSF) <r(s,u): 



g\sL) 



(2) 



It represents an integrated beta function for fi- 
nite scale transformations with rescaling factor s 
||. In principle, if one has control over the exact 
SSF, one can trace the non-perturbative evolu- 
tion of the coupling in discrete steps ~g 2 {L) — > 
g 2 (sL) — > g 2 (s 2 L) — > • • ■ . For one rescaling step, 
a is computed as the continuum limit a — > of 
the SSF at finite resolution, £(s,u, a/L). So far 
the method of construction for each u-value has 
been to choose L/a, determine (3 to match g 2 = u 
on the corresponding lattice and to then simulate 
on an sL/a-lattice with the same (3 and read off 
S = ~g 2 (sL). For further details we refer to 
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Fig. 1. Joint continuum extrapolation of E and 
E, and result of extrapolating E alone (triangle). 

3. SIMULATIONS AND RESULTS 

3.1. New simulation strategy 

In our previous computations L max defined by 
7? 2 (L max ) = 3.48 has been both an intermedi- 
ate reference scale and the low energy bound of 
our scaling investigations. In ref. || this scale 
was connected with the scale r$ « 0.5 fm from 
the static potential yielding L max /r — 0.718(16) 
in the continuum limit. In the same reference 
a relation connecting (3 with r /a was given for 
the range 5.7 < (3 < 6.57. These combined in- 
formations enable us to directly determine j3- 
values corresponding to L = s£ max for a series of 
chosen resolutions L/a for continuum extrapola- 
tion. We just use the cited relation for ro/a = 
(ro/L max )(L/a)(l/s) in the allowed range. In 
this way we construct an alternative lattice es- 
timator E(s, 3.48, a/L) with the same continuum 
limit <t(s, 3.48) but cutoff effects different from E. 
in fig. ry we demonstrate consistency between the 
two approaches for s = 1.5. In fig. |^ continuum 
extrapolations of E are shown for other s-values, 
which now take us down to 2.5L max in energy. 

For the 0(a) boundary improvement term of 
the SF we use the 2-loop approximation ||. For 
our Monte Carlo simulations we employed a sim- 
ilar 'hybrid-overrelaxation' algorithm as in |l]]. 
Each iteration consists of 1 heatbath update and 
iVoR subsequent overrelaxation sweeps (typically 
-^or = 5,7). As detailed there, non-gaussian 
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Fig. 2. Extrapolations of E(s, 3. 48, a/L), dis- 
carding the two coarsest lattices. 

tails and long autocorrelations in the coupling 
at low energy are overcome by a modified sam- 
pling procedure that enhances the tail contribu- 
tions and is compensated by rewcighting. 

3.2. Results 

The SSF E-values in fig. || approach the con- 
tinuum at rate compatible with a 2 . This shows 
that lattice artifacts linear in a, which in the SF a 
priori exist, are invisible to our accuracy with 2- 
loop boundary improvement and allowed us to fit 
by the form E = a + p 2 (a 2 /?"o) ■ As a safeguard 
against underestimating the uncertainty in a we 
omitted the two coarsest lattices. This is in con- 
trast with the extrapolation of L max /ro, where a 
small linear component was found pM . 
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Fig. 3. Running of a = o?sf = ff 2 /47r. 
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Fig. 4. L-dependence of ~g 2 versus PT. The 
artificial "4-loop" j3-function uses = l/(47r) . 
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Fig. 5. Estimates of m via leading-order strong 
coupling (full line) and our data (dotted line). 

Figs. H and |J display the scale evolution of the 
SF coupling over the whole energy range that 
is available from refs. |l]|| together with this 
work and confront it with results from integra- 
tions with the perturbative /3-function. With de- 
creasing energy a substantial deviation from PT 
now becomes clearly visible. The steep growth 
of ln{g 2 (L)} with L reveals the non-perturbative 
behaviour to set in for L > 0.7 ro. 

3.3. Discussion of the large —L asymptotics 

In the low-energy domain ~g 2 is dominated by 
non-perturbative contributions. Since the bound- 
ary fields are locally pure gauge configurations, 
any dependence of the effective action T on the 
background field is caused by correlations around 
the spatially periodic lattice. At large volume 
they are exponentially suppressed in a theory 



with a mass gap, and one expects 

g 2 {L) oc cxp{mL} at large L . 

The mass m is characteristic for the gauge field 
dynamics close to the boundaries (but not ob- 
viously related to the bulk correlation length). 
Based on this argument we attempt to compare 
the behaviour of ~g 2 as a function of L with the 
leading-order prediction from the strong-coupling 
expansion (llj: g~ 2 (L) ~ A x cxp{— m(/3)L}, 
with to = (3/4) tog = —3 In/? and toq the ++ 
glueball mass. Thus, assuming this ansatz, we fit: 

\n{g 2 (L)} =A'+ (m£ max ) x (L/L max ) . (4) 

Fig. H confronts the outcome of this analysis ap- 
plied to the two highest points with the slope ob- 
tained in leading-order strong coupling, 

TOr = (3/4) x (TO G r ) mL max w 2.3 , (5) 

with to-g^o w 4.3 from [12 ]. Exploring the pos- 
sibility, whether in eq. (H) a prefactor L c , c ^ 0, 
can build up when higher orders in the strong- 
coupling expansion are summed, would demand 
a calculation of at least a few further orders. 
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